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<D 1 ABSTRACT 

m ; _ . . . 

. We present a calculation of the rare decay modes D° — ► e + e 7 and D° — > 7 in the 

I framework of Standard Model. For the short distance part, we have derived QCD corrections to 

the Wilson coefficients involved, including Cg. The latter is found to be strongly suppressed by the 
corrections, leading to diminished values for the c — > ul + l~ branching ratios in the 10 -10 range. 
Within SM the exclusive decays are dominated by long distance effects. Nonresonant contributions 
are estimated using heavy quark and chiral symmetries to be at the level of 10%, compared to the 
contributions arising from D° — > V7 — ► l + l~j, with V = p,uj,4>. The total SM branching ratio is 



predicted to be in the range (1 — 3) x 1CP 9 . We also consider contributions coming from MSSM 
with and without R parity conservation. Effects from MSSM are significant only for the i?-parity 
violating case. Such contributions enhance the branching ratio D° — > ^ + /i~7 to < 0.5 x 10~ 7 , based 
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on appropriately allowed values for Cg and Cio- This selects D —>■ fi + fi 7 as a possible probe of 
new physics. 

PACS numbers: 13.25.Ft, 13.20.-v, 13.60.-r, 12.60.Jv, 



1 Introduction 



The charm physics is entering an exciting era. The high statistics and an excellent quality of data 
at FOCUS experiment now allow, among others, for high precision studies of charm semileptonic 
decays [1], determination of D 0,ziz decay times below 1% error level [2], as well as for searches of CP 
violation and rare D decays [3,4]. There is a very rich potential for charm physics at B-factories, 
with both Belle and BaBar having an active program in charm studies [5,6]. For instance, more 
than 120 milion charm pairs have already been produced at BaBar. This corresponds to more than 
220 000 D*-tagged D° decays, which will allow for precision lifetime and D° mixing analyses as well 
as for searches of rare charm decays [6]. An exciting charm physics program is under way also at 
CLEO, that was recently able to measure T(D*) for the first time [7, 8]. Among the rare D decays, 
the decays D -> V-/ and D -> V(P)l + l~ are subjects of CLEO and FERMILAB searches [9, 10]. 
In the following years a great phenomenological impact is expected from proposed CLEO-c physics 
programme. Next year more than 6 million tagged D decays are expected to be measured. This will 
allow for precision charm branching ratio measurements and consequently improved measurements 
of CKM matrix elements also in 6-sector, as well as for extensive studies of .D-mixing, CP violation 
and rare decays in the charm sector [11]. 

Parallel to the experimental studies, there has been an ongoing theoretical effort to understand 
charm physics. A number of studies has focused on the rare charm decays [12-20] and a possible 
impact of new physics on the predicted branching ratios [21-26]. Note however, that in rare D 
decays the nonperturbative physics of light quarks is expected to dominate the decay rates. Consider 
for instance the case of c — > wy transition that occurs only at one loop level in the Standard Model. 
The contributions coming from b, s, d quarks running in the loop are 

{0(X 5 mD : b — quark, 
0(\m 2 s ) : s- quark, (1) 
0(XAq CD ) : d— quark, 

where we have tentatively set Aqcd instead of m u for the u— quark contribution, anticipating 
the size of nonperturbative effects. The situation is quite different from the s — ► d FCNC (e.g. 
s — »■ dvv), where the same CKM hierarchy is present, but with the top quark replacing the b quark. 
Since b quark is much lighter than the top quark, it cannot surpass the A 4 suppression. Thus the 
contributions from the heaviest, 6-quark, are expected to be the least important. One can then 
expect that in rare D decays the nonperturbative long distance (LD) effects coming from the lighter 
two down quarks, d, s will give the dominant contributions. 

Since LD effects are difficult to control theoretically one would like to either find decay modes 
where LD effects are as small as possible and/or find observables where LD effects cancel. Such 
an observable was constructed in [23], where D° — > (p,a>)7 decays were considered. It was found 
that most of the LD effects cancel in the difference of the appropriately renormalized decay widths, 
making these channels a useful probe of new physics. Another interesting analysis is connected with 
decay modes D — ► (P, V)l + l~, where P = n,K, rj are the light pseudoscalar mesons, V = p,u, . . . 
are the light vector mesons and / + /~ is an electron or muon lepton pair. The decays have been 
estimated both in the SM and MSSM [17, 25-28]. In Ref. [25] it was found that the experimental 
bounds on Br(D + — > 7r + p + p~), Br(D° — > p fi + fi~) constrain the sizes of relevant trilinear R parity 
violating coupling more stringently than analyses from other processes. Measurements of rare D 
meson decays can thus already now constrain new physics scenarios in the up-like quark sector. 

In this paper we investigate the rare decays D° — > e + e~7, D° — ► [i + both in the Standard 
Model and in MSSM. A Standard Model analysis of D° — > branching ratios neglecting 
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QCD effects and LD transitions, has been made in Ref. [20], giving Br(D° — > = 6.3 x 

1CT 11 . However, LD effects are expected to dominate the SM prediction similarly to the D — > 
(P, decays. To evaluate the nonresonant LD effects we use the heavy quark effective theory 

combined with chiral perturbation theory (HQxPT) [29] . This approach was used before for treating 
D* strong and electromagnetic decays [30,31], as well as the leptonic and semileptonic decays of 
D meson (see [30] and references therein) and D° — > 77 decay [18]. In addition, we include 
contributions of vector resonances in our analysis. 

Another expectation based on the experience from D — > (P, V)l + l~ decays is that there are 
possibly large contributions in D° — > decays coming from SM extensions such as MSSM 

with R-parity violation. These expectations make D° — > channels interesting from both 

experimental as well as from theoretical side. 

Our calculations show that as a result of the LD contributions, these decays would occur with 
a branching ratio of (1 — 3) x 10~ 9 in the SM, nearly two orders of magnitude larger than in the 
previous estimate [20]. Moreover, MSSM with P-parity violation as presently restricted, allows for 
a branching ratio of D° — > in the 10 -7 range. 

The paper is organized as follows. We start with the Standard Model prediction in section 2, 
where first a discussion of renormalization group improved effective weak Lagrangian together with 
the calculation of c — ► ul + l~ inclusive mode is given. This is then followed by the estimates of 
nonresonant as well as of resonant LD contributions to the decay width D° — > In section 3 

we present possible effects of SUSY extensions of Standard Model. In appendices we collect some 
further details about the calculation of c — > u effective weak Lagrangian, as well as the explicit 
formulae of the calculations. 



2 Standard Model calculation 

We will devote the first part of the present paper to the estimation of D° — ► decay width in 

the context of the Standard Model. At the quark level, this decay mode cannot proceed through 
tree diagrams and is thus induced only at the one loop level in the Standard Model. Possible quark 
diagrams are shown on Fig. 1. These then translate into an effective weak Lagrangian at the scale 
of m c . 



2.1 Effective weak Lagrangian 

The effective Lagrangian describing the weak c — > u transitions at the scale of \i = m c is (see 
Appendix A, Eqs. (A.4)-(A.6)) 

C eS = -^[v: d V ud J2CiQf + V c * s V us J2c i Qi-V^V ub C iQi\> ( 2 ) 

* i=l,2 i=l,2 j=3,...,10 
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where 



Ql = (u a qP) v „ A (fc a ) v „ A , 
Q3 = (uc) V - A ^2(qq)v-A, 



Q\ =(uq)v-A(qc)v-A, 



Q 4 =(u a c^) V - A J2(fQ a )v-A, 



(3a) 
(3b) 



Q5 = (uc) V - A ^2(qq)v+A, 



Q 6 =(« a c\_ A ^(/ g a ) m , 



(3c) 




Qi = -^mcF^ua^ P R c, 
e 2 

Q9 = --^{ul^cl)^^), 



(3d) 



(3e) 



with qL = Piq and Pr,l = (1 =t 7s)/2 the chirality projection operators, while we have suppressed 
the color indices in the currents of form (qq 1 ) = (q a q' a ). The sum over q runs over the active 
quark-flavors. At scale fi ~ m c these are Wilson coefficients to which QCD 

corrections are administered. We do not include in the analysis higher dimension operators. 

Note, that the penguin operators Q3 10 are proportional to V* b V u b matrix elements in the 
effective weak Lagrangian (2). In Wolfenstein parametrization this is ~ A , which has to be 
compared to the CKM suppression of Qi$ operators, V CS V US ~ A, where A = sin6* c = 0.22. Penguin 
operators are thus greatly suppressed in AC = 1 transitions. They are relevant only in special 
observables such as CP asymmetries [32]. In the literature [20,25,26] as an estimate for Cg(/j, c ) 
Wilson coefficient, the result from electroweak theory without QCD, Cg L (where IL stands for 
Inami-Lim [33]) has been used. Since Cg L is not V u i, suppressed, it greatly overestimates the effect 
of Qg operator insertion on the predicted decay widths. We will thus devote the rest of this section 
to clarify this point. 

The values of Wilson coefficients C±, . . . , C10 at scale fx = m c are obtained by using the same 
method as in the existing calculations for s — ► d transitions [34-36] at leading (LO) and next-to- 
leading order (NLO). Application to c — ► u transition is straightforward, but some care has to be 
taken when integrating out the 6-quark at the intermediate step of renormalization group (RG) 
evolution. The charge of the intermediate 6-quark is important for the matching of electroweak 
Cg Wilson coefficient. Since this calculation has not yet been performed we give further details 
in appendix A. The Wilson coefficients C\, . . . ,Cq for c — > u transitions have been calculated 
already in [32] at NLO, while the LO calculation of C\ has been presented in [12, 14]. The values 
of calculated Wilson coefficients are listed in Table 1. For a comparison the values of Wilson 
coefficients at LO order are given as well, but calculated with the two-loop evolution of the strong 
coupling constant. The values are given for the central value of A^ 5 ) = 216 ± 25 MeV and the 
matching scale = 4.25 GeV. The one sigma change in A^ 5 ) corresponds to a change of about 
10% in Ci v .. 5 6. We find a pronounced scale dependence for the Cg coefficient below 1.5 GeV, as a 
consequence of large cancelations in the RG evolution equations. The situation is very similar to 
the case of coefficient Z-jy in Kl — > 7r°e + e _ [34]. The LO value of Cg even changes sign near fi ~ 1 
GeV, being positive for [i > 1 GeV. Note, that uncertainties in the value of the Cg coefficient will 
not propagate into the decay rates as the Qg operator is V u b suppressed. Note also, that Qio does 
not mix with other operators due to chirality, so that Cio(/i c ) = Cio(Miy) ~ 0. 

As for the C7 Wilson coefficient, the leading order mixing of operators Qj t % with operators Qi,...,6 
vanishes. It is only at two-loop that the anomalous dimension matrix has nonzero values mixing 
Ci,...,6 into C7. Since two-loop results are scheme dependent, it is then customary to introduce 
an effective anomalous dimension matrix 7(°) cff [36], which is scheme independent as is the case 
in leading order results. Using LO anomalous dimension matrix j^ cS and NLO evolution for a s , 
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Ai(GeV) 


Ci 


c 2 


c 3 


c 4 


c 5 


c 6 


c 9 


LO 


1.0 


-0.64 


1.34 


0.016 


-0.036 


0.010 


-0.046 


-0.07 


NLO 


1.0 


-0.49 


1.26 


0.024 


-0.060 


0.015 


-0.060 


-0.60 


NLO 


1.5 


-0.37 


1.18 


0.013 


-0.036 


0.012 


-0.033 


-0.13 


NLO 


2.0 


-0.30 


1.14 


0.009 


-0.025 


0.009 


-0.021 


-0.13 



Table 1: Values of Wilson coefficients at scales n = 1.0 GeV, 1.5 GeV, 2.0 GeV, calculated at next-to-leading order 
(NLO) as explained in text. For a comparison in the first line the LO values are given at scale \x = 1 GeV, but 
calculated with two loop evolution of strong coupling constant. 



m b = 4.25 GeV, we arrive at (see also [14]) 

Cf (1.0 GeV) = 0.13, Cf (1.5 GeV) = 0.087, Cf (2.0 GeV) = 0.065, (4) 

Note that, as we already mentioned, the Wilson coefficient Cg(fi c ) has been estimated previously 
[20,25,26] by using the result from electroweak theory without QCD, i.e. taking Cg L , based on 
the (unproved) expectation that Cg is not much affected by QCD corrections. The leading order 
expression in terms of 'm ds /rn'^ 1S 1 

Cf ~ - A s 16/9 In (m s /m d ), (5) 

where Xj = V^Vuj / (V* b V ub ) . Using m s /m d = 17-22 [37] we arrive at the value V* b V ub Cf^ ~ 
-U c *K s 16/91n (m s /m d ) = -1.13 ± 0.06 which should be compared to V* b V ub Cg(n) ~ 10~ 4 . The 
value of Wilson coefficient is thus four orders of magnitude smaller than the corresponding parame- 
ter obtained by neglecting QCD interactions! The reason for this discrepancy lies in the appearance 
of large logarithms \n(m^ s /mw) that avoid the GIM suppression otherwise present in Cg. It is 
exactly these large logarithms that RG evolution sums correctly. Since small scales of order m dtS 
lie in the nonperturbative region of QCD, the approximation of using (5) without QCD corrections 
is not valid. 

The logarithm appearing in (5) is exactly reproduced in the calculation of inclusive modes 
c — ► ul + l~ , if mass-independent renormalization is used (see appendix C of [38]). To show this 
explicitly, we consider the calculation of c — ► ul + l~ in the naive dimensional regularization (NDR). 
The amplitude can be parametrized as 



Gf , 
y/2 



M = -^=v: b v ub 



Cf(Q 7 }° + Cf(Q 9 )° + CT ff <Q 10 > l , (6) 



with (<57,9,io)° the tree level matrix elements of the operators. Note that 10 are not Wilson 
coefficients but merely parametrize the invariant amplitude. The Cg S coefficient is dominated by 
the 1-loop contributions coming from insertion of Q\ 2 operators, q = d,s. The virtual photon is 
emitted from the intermediate d, s quarks. This contribution is of order a® and proportional to 
V* q V uq and is thus only once Cabibbo suppressed. Using existing results for b — > sl + l~ at NLO 
[39-41], we arrive at 

V; b V ub Cf = 2V: S V US (h(z s , s) - h(z d , J)) (3d(m c ) + C 2 (m c )) , (7) 



x For further details about the calculation see appendix A, where also a discussion regarding 6*7,10 is presented 
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with z q = m q /m c , s = (m;+;-/m c ) 2 and mi+i- the mass of the lepton pair, while 



, / v 8 8 4 

H z > 8 ) =~ g lnz+ 27 + g x 



in 



-(2 + x) v / |l-x| <( 



2 Arctan 



27T, 



1 



for x < 1, 
for x > 1, 



(8) 



where x = 4z 2 /s. In (7) the contributions suppressed by V* b V u b are neglected. These include 
the tree-level contribution from Qg as well as 1-loop contributions coming from insertions of QCD 
penguin operators Q3,...6- From expression (7) one should reproduce the Inami-Lim result (5), when 
momenta and masses of external particles are set to zero. Taking the limit m^- <C md :S , one gets 



lim (h(z s , s) 



h(z d ,s)) 



- In (— \ 

~9 n \m d J 



(9) 



Taking the values of C\p, Wilson coefficients at the weak scale C\ ~ 0, C2 ~ 1 one arrives at 
the Inami-Lim result (5), as expected. Note, that the logarithm ln(m^/m s ) in (5) arises from the 
insertion of Qi,2 operators. Phenomenologically more interesting is the limit mj+j- ~ m c 3> niac- 
in the limit mj+j- — ► 00 the difference (h(z s , s) — h{z d , s)) vanishes, while for mi+i- ~ m c it is at a 
level of few percent! Using [25, 26] Cg L (5) instead of Cg S (7), which includes the QCD corrections, 
one overestimates the dBr(c — ► l + l~)/ds. 
Explicitly, the branching ratio is [26] 



Br(c -► «/+/-) 



"768vr 5 r(£> ) 
+ 



|y c * 6 K 6 | 2 (i-5) 2 



1 + Ts [\cf\ 2 + \C$\ 2 ) + (c 7 en *C 9 



1 + 

cff* AcfT 



cff|2 



16 



(10) 



where we write s = (mi+i- /m c ) 2 as before. For the value of Cj S we use the two-loop result of Ref. 
[14], Cf = A s (0.007 + 0.020i)(l ± 0.2), with X s defined after Eq. (5). The dominant contribution 
to C| ff comes from the insertion of Q\ operator, while the contributions from the insertion of Q\ 
operators vanish because of color structure. The coefficient Cf^ — in the Standard Model. 
Using m c = 1.4 GeV one arrives at 



Br(c — ► ue + e 
Br(c — ► ufi + fi~ 



2.4 x 10~ 10 , 
0.5 x 10~ 10 , 



(11) 



where the dominant contribution comes from the part of the amplitude. This is in contrast 
to Refs. [25,26], where Cf was estimated using Cg L . This lead to the branching ratios of one 
(for e + e~) to two (for fi + orders of magnitude higher, with Cg contribution dominating the 
branching ratio. 

The suppression of QCD corrected Cf (7) compared to C£ L (5) comes from two sources. The 
cancellation of s and d quark contributions in (7) is very strong even at moderate values of s, with 
(h(z s ,s) — h(zd,s)) < 10% for s > 0.3. There is also a sizable cancelation between C\(m c ) and 
C2(m c ) in (7). This cancelations could in principle be modified by the two- loop QCD corrections 
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SD.l 



SD.4 




Figure 2: The short distance diagrams. The effective weak Lagrangian vertex is denoted by an empty square. The 
relevant operator is denoted as well. 



to Qi,2 matrix elements 2 . If the cancelations were completely lifted, one can estimate the possible 
effect by Cf {m c )Cl L . This leads to roughly the same prediction for Br(c — > ue + e ), while it 
can increase Br(c — > u/i + /[/~), as Cg ff affects mostly the higher s part of the decay width distribution. 

Note, that the calculation of c — > ul + l~ is in many respects different than the calculation of 
6 — ► sl + l~ . The operators Q"'^ 6- ^ in 6 — > are equivalent to Q12 operators in the c — > 

transition, but with different CKM factors multiplying the operators in the effective Lagrangian. 
In b — > sl + l~ then only the Q^^p*^ operators contribute, as the contributions coming from the 

Qi*2~~* S ^ operators are V^j, suppressed. Hence, there is no approximate cancellation of the type 
(h(z s ,s) — h(zd,s)) found above. Note also, that in b — > sl + l~ the penguin operators Q3,...,io are 
not CKM suppressed relative to Qi,2 and have to be taken into account, contrary to the c — > 
case, where penguin operators are suppressed. 

The V u b suppression of penguin operators Q3 10 is of course present in the calculation of 
exclusive charm decays, where the insertions of Q\p operators again dominate the rate. This will 
be discussed in more detail for the case of D° — ► decay in the following section. Before 

we proceed with the calculation, let us mention the commonly used terminology of long distance 
(LD) and short distance (SD) contributions. These are usually separated in the discussion of weak 
radiative decays q' — > (777 or q' — > (77 decays. The SD contribution in these transitions is a result 
of the penguin-like transition induced by the operators (37,9,10, while the long distance contribution 
arises from insertions of Qi,2 operators, when the off- or on-shell photon is emitted from the quark 
legs. We will follow this classification in the following. 

At this point, we mention our result for the SD contribution coming from the operators (57,9,io 
(see Fig. 2). This contribution turns out indeed to be very small in the SM, due to the CKM 
suppression (2). Evaluating the expectation values of operators ^7,9,10 by using heavy quark 
symmetry as described in Eq. (23) of the next sub-section (see the explicit expressions in Appendix 
D) and using the values of Wilson coefficients listed in Table 1 and in Eq. (4), one arrives at the 
corresponding branching ratios for D° — ► / + / _ 7 of 1CT 17 — 1CT 18 . This is negligible compared to 
the LD contributions calculated in the subsections 2.2, 2.3. Our result for the SD contribution to 
these decays is several orders of magnitude smaller than the result of [20] , which was obtained with 
an unrealistic value of Cg that did not include QCD corrections. 

2 The existing two-loop calculations of Qi,2 matrix elements in b — > sl + l~ [42,43] have been done for small s, 
where no substantial increase in c — » ul + l~ is expected. 
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2.2 Nonresonant LD contributions 



Turning now to the LD contributions, we start with the nonresonant contributions. As we will see 
later on, it is in the nonresonant contributions that the extensions of Standard Model can show up. 

The most general invariant amplitude for D" — > Z + Z~7 decay following from the effective La- 
grangian (2) is 3 



M = M^e;(k)^u( Pl ) lu v( P2 ) + M^el{k)^u{jp 1 ) lv ^v{jp 2 )+ 



+ M BS (p 2 )[u( Pl )( 



f j>D _ p'd f 
. pi ■ k p 2 ■ k 



~)l5v{p 2 )], 



where 



Kl = Co, 5 (p 2 )(v» u - ^r) + D , 5 (p 2 )e^k aP(3 , 



(12) 



(13) 



with pi£ the four-momenta of lepton and antilepton respectively, p = p\ + p 2 the momentum of 
lepton pair, k the photon momentum and its polarization vector. The form factors Co^(p 2 ), 
Dq^(p 2 ), Mbs{p 2 ) are functions of p 2 only and in particular do not depend on k-p\ or k-p 2 . Co,D§ 
are parity violating terms, while C^,Dq and the bremsstrahlung part of the amplitude, Mbs, are 
parity conserving. 

The partial decay width is then 



k-p 



dT _ 1 
dp 2 lQ7r 3 m 3 D 1 3p 



(|C | 2 + |Z?o| 2 (M 2 )(l + 2Ap) + 



+ (|C5| 2 + | j d 5 | 2 (M 2 )(i-4A') 



+ ^ [ ((P 2 ) 2 + m 2 D (m 2 D - 4m 2 )) In (— 



2p 2 m 2 D 



(14) 



m 



+ 4Z(D M^ s )-(k- P y In 



p* 



where ft 2 = m 2 /p 2 , with m the lepton mass, ^ = ^/l — 4fi 2 , while k-p = {m D —p )/2. We checked 

that this expression agrees with the similar expression for the partial decay width Kl — »■ as 
given in [44, 45], as well as with the B — > decay width as given in [46]. 

The nonresonant LD contributions will arise in our approach from chiral loop contributions 
shown on Figure 3. The weak vertices receive contributions from Q\ t2 operators in the effective 
Lagrangian (2). The sizes of these contributions are estimated using factorization approximation. 
The effective [47] four quark nonleptonic AC = 1 weak Lagrangian is then 



G F 



~m E V W jV^[a 1 (^c)(ur^q) + a 2 (ur' i c)(qr il q) 
v 2 „f-f„ 



(15) 



q=d,s 



where V 1 = 7 M (1 — 75), ai are effective Wilson coefficients, V qiqj the CKM matrix elements, while 
products of currents in (15) are understood to be evaluated in the factorization approximation. 
We use the phenomenologically motivated values 4 01 = 1.26, a 2 = —0.49 of "new factorization" 



3 We use e 0123 = 1. 

4 The new factorization values of effective Wilson coefficients correspond to iV c 
Wilson coefficients ai = C2, a 2 = Ci. 



00 limit and are in terms of 
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M 2J 



fef 4 _l ^X^ M 4 _2 M 4-3 ^X^ M^£~£ M 4 _g V-^{^ M 4 _ 6 

^ =8® ' ^ ® ® ' ® ^ ' ®® ' ® ® ' ® ® 



Figure 3: Nonvanishing one loop diagrams. The dashed lines represent charged Goldstone bosons flowing in the 
loop (K + , 7r + ), while the double lines represent heavy mesons, D and D* . The two crossed circles denote weak vertex 
calculated in the factorization approximation. The sum of diagrams in each row is gauge invariant and finite. 



[48]. It is worth pointing out that long distance interactions will contribute only if the SU(3) 
flavor symmetry is broken, i.e. if m s / m^. Namely, due to V U( iV* d ~ — V US V* S , if = m s the 
contributions arising from the weak Lagrangian (15) cancel. Note also that in diagrams of Fig. 3 
only the term proportional to a\ contributes. The 02 part of effective Lagrangian (15) gives rise to 
the resonant LD contributions and will be discussed in the next section. 

We calculate the nonresonant LD contributions in the framework of Heavy Quark Chiral Per- 
turbation Theory HQ%PT [18, 30]. This model will serve us when hadronizing the currents [29] of 
the quark effective weak Lagrangian. In the framework of HQxPT a number of coupling constants 
appear that are fixed from experiment as discussed in Ref. [18] and are listed in Table 2. In the 
following we first give a short introduction to HQ^PT and then turn to the discussion of results. 

In the leading order of HQxPT the light pseudoscalar mesons are described by the usual 0(p 2 ) 
chiral Lagrangian 

C%> = ^ tr^E^Et) + tp- tr(M,£ + M q tf) , (16) 

where E = exp(2ill//) with LT = Y2j -7^ A- 7 7r J containing the Goldstone bosons -ir,K,r), while the 
trace tr runs over flavor indices and AAq — diag(?77» u , md,m s ) is the current quark mass matrix. 
From this Lagrangian we can deduce the light weak current of the order 0(p) 

il = -^tr(S^StA a ), (17) 

corresponding to the quark current = qL7n^ a QL (with A a an SU(3) flavor matrix). 

For the heavy mesons interacting with light pseudoscalars we have the following lowest order 
0(p) chiral Lagrangian 

62 = - Tr(H a ivD ab H b ) + g Tr(H a H bl ^ a 75) , (18) 

where D^ b Hh = d^Ha — H b V ba , while the trace Tr runs over Dirac indices. Note that in (18) and 
the rest of this section a and b are flavor indices. The vector and axial vector fields V M and A^ in 
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(18) are given by: 

V M = \{td^ + gd£), Ap = \{£}d^ - id^) , (19) 

where £ = exp (ill//). The heavy meson field H a contains a spin zero P a and spin one P afl boson 
fields 

H a = P + (P^-P alb ), H a = 7 °(H a )h°, (20) 

with P± = (1± y)/2 the projection operators. 

From symmetry grounds, the heavy-light weak current is bosonized in the following way [29] 

q a ^P L Q = Y Tr[^ P L H b ^J, (21) 

where Pr^l = (1± 7s)/2, Q is the heavy quark field in the full theory, in our case a c-quark field, 
and q is the light quark field. Note that the current (21) is O(p ) in the chiral counting. The 
constant a is related to the physical decay constant fr> through the well known matrix element 

(0\u 7 ^ 5 c\D°) = ip%f D , (22) 

from which a = f D y/m~5- From [37] one deduces f Ds = 268 ± 25 GeV and a = 0.38 ± 0.04 GeV 3/2 . 
In the same way as the heavydight current (21), operators of more general structure (uTc), with T 
an arbitrary product of Dirac matrices, can be translated into an operator containing meson fields 
only [49] 

(uFc) - | TT[P R TH4 a ] + ^ Tr[P L TH b ^ ba }. (23) 



/ 


132 MeV 


9 


0.59 ± 0.08 


a 


0.38 ± 0.04 GeV 3/2 


P 


2.3 ± 0.2 GeV" 1 


ai 


1.26 


02 


-0.49 



Table 2: Coupling constants appearing in HQxPT, that is used in the estimates of nonresonant contributions. For 
further details see text and [18]. In the last row values of effective Wilson coefficients are given [48]. Loop integrals 
are calculated in MS scheme with scale n — 1 GeV, while in (24) m c = 1.4 GeV. 



The photon couplings are obtained by gauging the Lagrangians (16), (18) and the light current 
(17) with the U(l) photon field B^. The covariant derivatives are then T>^ b H b = d^H a +ieB^(Q' H— 
HQ) a - H b Vg a and = + ieB^Q,^] with Q = diag(§,-±, -±) and Q' = § (for the case 
of c quark). The vector and axial vector fields (19) and the light weak current (17) contain after 
gauging the covariant derivative T>^ instead of d^. However, the gauging procedure alone does not 
introduce a transition DD*^y without emission or absorption of additional Goldstone boson. To 
describe this electromagnetic interaction we follow [50] introducing an additional gauge invariant 
contact term with coupling of dimension -1 

Cf3 = -^Ti H a H b a^F^Ql a - -£-Q' Tr H a a^H a F^, (24) 

where Q* = §(£ f Q£ + £Q&) and = d^B v - d v B^. The first term concerns the contribution of 
the light quarks in the heavy meson and the second term describes emission of a photon from the 
heavy quark. Its coefficient is fixed by heavy quark symmetry. From this "anomalous" interaction, 



9 



Figure 4: One particle reducible diagrams with photon connecting initial (pseudo)scalar and final state particles 
are zero. 



both H*H~f and H*H*j interaction terms arise. Even though the Lagrangian (24) is formally 
1/rriQ ~ m q suppressed, we do not neglect it completely. We do not take it into account in chiral 
loop contributions of Fig. 3, as it has been found to give a rather small contribution in a very similar 
case of D° — ► 77 analysis [18]. The D° — ► D°*j transition will be, however, needed to estimate the 
short distance contributions shown on Fig 2. These will give numerically irrelevant contributions 
for SM predictions but will be important later on, when we extend the analysis to MSSM case. 
Note also, that the Lagrangian (24) in principle receives a number of other contributions at the 
order of 1/mg. However, these can be absorbed in the definition of (3 for the processes considered 
[50]. 

Using HQ^PT as described above, one arrives at the set of nonzero 0(p 3 ) dia grams listed in 
Fig. 3. Each row of diagrams on Fig 3 is a gauge invariant set. The sum of diagrams in each 
row is also finite. Separate diagrams are in general divergent and are regulated using dimensional 
regularization. Further details on this subject can be found in the appendix B. The explicit 
expressions of the corresponding amplitudes can be found in appendix C. Note that the chiral loop 
contributions of Fig. 3 contribute only to the Mq V part of the invariant amplitude (12) . Namely, 
the l + l~ pair couples to the charged mesons in the loop only via electromagnetic current. This 
also leads to the 1/p 2 photon pole in the amplitude (p being the momentum of the lepton pair). 
The LD nonresonant contributions coming from Fig. 3 thus exhibit a pole behaviour at small 
lepton momenta. This pole is either cut off by the phase space because of nonzero lepton masses 
(p 2 = 4m 2 ), or by experimental limitations due to Dalitz conversion [25]. 

Note that there is no photon bremsstrahlung off the final lepton pair in the chiral loop contri- 
butions. Namely, diagrams of the type shown on Fig. 4, with initial meson being a (pseudo) scalar, 
and with a photon connecting the two blobs, vanish due to gauge invariance. 

The diagrams of Fig. 3 are evaluated in the minimal subtraction (MS) renormalization scheme. 
However, the sum of diagrams is finite and scheme independent. We use the values of coupling 
constants listed in Table 2. Integrating over the whole available phase space one arrives at the 
estimates 

Bl (D° -► e + e-7) nonrcs = 1.29 x 1(T 10 , Br(D° -► fi + M _ 7)nonr es = 0.21 x 1(T 10 . (25) 

Due to a photon pole, the larger part of the electron channel branching ratio comes from the region 
of the phase space with p 2 ~ 0. The phase space is cut off by muon masses at much higher p 2 , 
giving a smaller contribution of nonresonant LD effects to this decay channel. 

2.3 Resonant LD contributions 

The mechanism of the decay D° — ► through resonant intermediate state is depicted on Fig 5. 

The D° meson first decays into a vector meson and a photon, D° — > Vj. The vector meson than 
decays into a lepton pair, completing the cascade D° — > Vj — > The decay width coming 
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Figure 5: The mechanism of D° — ► l + l 7 decay through intermediate vector resonance state V. 

from this mechanism can be written as [51] 

dV D o^ Vl ^i+i- 1 1 - 



dp 2 



■ V "<Tr(M? 



V 



p 2 ) 2 + M 2 T 2 



(26) 



where p is the momentum of the lepton pair, while My and V are the mass and the decay width of 
the vector meson resonance. Several assumptions go into the derivation of the simple, but physically- 
well motivated formula (26). First of all the interference with other channels is neglected. Under 
this approximation the formula is generally valid for the case of scalar resonances. Following the 
reasoning of Ref. [51] it is easy to show, that Eq. (26) is valid also for the case of electromagnetic 
decay of vector resonance into a lepton pair. 

Since vector resonances p,co,(fr are relatively narrow Eq. (26) can be further simplified using 
the narrow width approximation T <C My 



Br(D° -». U 7 Z+r 7 ) = Br^ -► V 7 )Br(U -► 



(27) 



The narrow width approximation is valid at 5% level for p, and below 1% for u, <f) mesons. To 
obtain numerical estimates, the experimental data on the branching ratios Br(V — ► l + l~) [37] can 
be used. On the other hand none of the decays D° — ► U 7 have been measured yet. We thus use the 
theoretical predictions of branching ratios Br(D° — > Vj). As the central values we use the recent 
predictions of Ref. [52], where a reanalysis of Ref. [16] has been performed using the quark model 
to determine relative phase uncertainties. As a comparison we also list in Table 4 the predictions 
of Ref. [12]. Note that for the upper limit predictions in [12] VMD model was used, with the 
main numerical input the experimental value of Br(D° — > p°4>). However, the central value of this 
branching fraction as cited in [37] has decreased by a factor of three between 1994-2002. Thus the 
upper limits on predictions of [12] should be divided by three, bringing the values in fair agreement 
with [52]. 



Decay 


Exp. [37] 


Decay 


Exp. [37] 


Br(p u e + e~) 


(4.54 ±0.10) x 10~ 5 


Br(p u -► p+p-) 


(4.60 ± 0.28) x 10" 5 


Br(oj — ► e + e _ ) 


(6.95 ±0.15) x 10~ 5 


Br(w — > p + p~ ) 


(9.0 ±3.1) x 10" 5 


Br(0 — ► e + e _ ) 


(2.96 ± 0.04) x 10" 4 


Br^^p+p-) 


(2.87^:S) >< 10-4 



Table 3: Branching ratios of vector mesons decaying to a lepton pair as compiled in Ref. [37]. 



Using the values compiled in Tables 3, 4 together with Eq. (27) one immediately arrives at 

Br(D° -> py -> ~ 5 x 10~ n , (28) 



Br( J D° 
Br(D° 



4>i 



/+r 7 ) 



x 10 



-11 



10" 



(29) 
(30) 
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Decay 


Theor. [52] 


Theor. [12] 


Exp. [37] 


Br(D u -> /) u 7) 


1.2 x lO -8 


(1 - 5) x 10 -0 


< 2.4 x 10" 4 


Br(L> u -► ury) 


1.2 x i(r B 


~ 2 x 10 -B 


< 2.4 x 1(T 4 


Br(L> u -> 7 ) 


3.3 x 10-° 


(1 - 34) x 10~ e 


< 1.9 x i(r 4 



Table 4: Theoretical predictions for decays D° — » Vj [12,52]. Predictions of Ref. [52] are used as central values 
(see also comments in text). In the last column the experimental upper limits are listed. 




0. 0.25 0.5 0.75 1. 1.25 1.5 1.75 0. 0.25 0.5 0.75 1. 1.25 1.5 1.75 

m ee [GeV] m„ M [GeV] 



Figure 6: The normalized decay width distribution (dF / dp 2 ) /T as function of effective lepton pair mass m l + l - 
(where wif+j- = p' 2 ) for e + e~ (left plot) and fi + pT (right plot) final lepton pair. The dotted line denotes SM 
nonresonant contribution, solid black line denotes full SM prediction, while dashed line denotes largest possible 
MSSM contribution with R parity violation. 



with l + l~ = e + e~ , fi + fi~ . Above we have used the fact that differences between e + e~ and /i + 
decay modes in the standard model come from the phase space differences only. These are relatively 
small compared to other theoretical and experimental uncertainties entering predictions (28)-(30), 
and are as such neglected. 

As seen from the estimates (28)-(30) the largest contribution to D° — ► Z + Z _ 7 comes from the 
intermediate eft resonance, being approximately one order of magnitude larger than the other two 
contributions. Note also, that in the region of p 2 , where vector resonances are important, the 
nonresonant contribution calculated in the previous section is several orders of magnitude smaller. 
We can thus safely neglect possible interference between nonresonant and resonant contributions 
and simply add resonant contributions (28)-(30) to the nonresonant ones (25). The decay width 
distribution is plotted on Fig. 6, while the predicted branching ratios are 

Bi(D° -» e + e -7) SM = 1.2 x 1(T 9 , Br(D° -» m + m~7)sm = 1.1 x 10~ 9 . (31) 

Note that if the values of Ref. [25] had been used, the predicted branching ratios could be utmost 
a factor of three higher. 

Incidentally Fig. 6 also explains, why the D° — > — > 7*7 cascade could be neglected in 
the D° — ► 77 decay rate calculation of Ref. [18]. Namely, for 7* almost on-shell the decay width 
is dominated by the nonresonant contributions. In the calculation of D — > 77 [18] these were 
described using HQxPT along the lines presented in section 2.2. 
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3 Beyond the Standard Model 



In this section we will consider possible effects of physics beyond the Standard Model, that could 
enhance the predicted branching ratios (31). The effects of new physics show up in the models we 
considered in the values of Wilson coefficients 



Cf™ = d + 5C U (32) 

where Cj are the SM values of Wilson coefficients listed in Table 1 and in Eq. (4), while 5Ci denote 
the changes due to new physics effects. Note that the general feature of all the SM extensions is 
to overcome the V* b V ub suppression of penguin operators (57,9,10 (2). Another general feature is 
that the new physics effects will extend the basis of penguin operator (3) by operators Q' 7 9 10 with 
quark chiralities switched (i.e. they are obtained by exchanging Pr <-> Pl in (3d),(3e)). 



3.1 Minimal Supersymmetric Standard Model 

We start with the simplest supersymmetric extension of SM, the Minimal Supersymmetric Standard 
Model (MSSM). It is constructed by putting the SM fermions in chiral multiplets and the SM gauge 
bosons in the vector multiplets, thus in effect doubling the spectrum of Standard Model fields. If 
no particular SUSY breaking mechanism is assumed the MSSM Lagrangian contains well over 100 
unknown parameters. It is thus very useful to adopt the so-called mass insertion approximation. 
In this approximation the basis of fermion and sfermion states is chosen such that all the couplings 
of these particles to neutral gauginos are flavor diagonal, but then the squark mass matrices are 
not diagonal. The squark propagators are then expanded in terms of nondiagonal elements, where 
mass insertions induce changes of squark flavor [53]. The mass insertions are parametrized as 

/ AM _ ( M ij) 2 AB ,„o\ 
q 

where i 7^ j are flavor indices, A, B denote chirality, (M-j) 2 are the off-diagonal elements of up-type 
squark mass matrices and Mq is the average squark mass. 

The largest contribution to c — > ul + l~ transition is expected from gluino-squark exchanges 
[25, 26, 54]. Allowing for only one insertion, the contributions from gluino-squark exchange diagrams 
are 



fJ{U , Pl32(z) Mg 
{0l2)LL "j h {0 12 )LrPi22{Z) 



(34a) 



V: b V ub 6C 9 = ^^^na^hLPo^z), (34b) 
V: b V ub SC w ~ 0, (34c) 



where z = M-/M-, while the functions Pijk(z) are 

Pijk(z) = [ dx f " - ^\ h - (35) 
Jo 



(1 — x + zx) k 



The Wilson coefficient C\q receives first nonzero contributions from double mass insertions, there- 
fore we neglect it in the following. The Wilson coefficients C 7 910 corresponding to the operators 
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with "wrong chirality" receive contributions from gluino-squark exchanges that are of the same 
form as expressions (34), but with the interchange L <-> R. 

In numerical evaluation of possible MSSM effects we use gluino mass M g = 250 GeV and the 
average squark mass Mq = 250 GeV, that are given by the lower experimental bounds [37]. For 
the bounds on the mass insertions we use the analysis of [24, 26]. The strongest bounds on mass 
insertion parameters (Sf 2 )LR are obtained by requiring that the minima of the scalar potential do 
not break charge or color, and that they are bounded from below [24, 55], giving 

KS12) lr\,\(S^)rl\ < 4.6 x 10- 3 , for M q = 250 GeV. (36) 

The bounds on mass insertions {5\ 2 )ll and {8\ 2 )rr can be obtained from the experimental upper 
bound on the mass difference in the neutral D system. Saturating the experimental bound Amp < 
4.5 x 10" 14 GeV [56, 57] by the gluino exchange gives [24, 26, 58] 

\{5i 2 )ll\ < 0.03, for M g = Mq = 250 GeV, (37) 

where (5i 2 )rr has been set to zero. These translate into 

\V* b V ub SC 7 \ < 0.04, \V* b V ub 5C' 7 \ < 0.04, (38a) 

\V: b V ub SC 9 \ < 0.0016, \V* b V ub 8C' 9 \ ~ 0. (38b) 

Note that both C7 and C' 7 receive largest contributions from {5\ 2 )lr insertions. Note also that the 
upper limits on C7 coefficient is three orders of magnitude larger than the Standard Model value, 
while for C9 is an order of magnitude larger than the SM value. However, as discussed in previous 
section, SM prediction is dominated by Q\^ 2 insertions and therefore by long distance effects. 

The contributing diagrams are shown on Fig. 2, to which the diagrams with Qi — > Q\ should 
be added. In mass insertion approximation the coefficient C\q is small and will be neglected in the 
following. When Q-j$ operators are inserted, the photon bremsstrahlung off the final lepton pair 
is not possible. In the case of Q7 operator this is because the diagrams are of the type shown in 
Fig. 4, while in the case of Qg operator the bremsstrahlung is prohibited because of vector current 
conservation. 

Taking the values of induced Wilson coefficients at the upper bounds we obtain 

Br(D° - e + e- 7 )MSSM = 1-4 x 10" 9 , Bi(D° - /xV~7)mssm = 1-2 x lO" 9 . (39) 

The MSSM contribution to the decay rate is entirely due to gluino exchange enhancement of C7, C' 7 
coefficients. The decay rate is thus enhanced in low p 2 region, which also explains larger increase 
of D° — ► e + e~7 decay rate. The increase is, however, not significant enough to dominate over the 
resonant contributions (28)-(30). MSSM effects, if any, are thus too small to be unambiguously 
detected experimentally in the decays D° — > 



3.2 R parity violation 

The situation is quite different once the assumption of R parity conservation is relaxed and the 
soft symmetry breaking terms are introduced. We follow the analysis of Ref. [25]. The tree level 
exchange of down squarks results in the effective interaction 

£eff = ^^^CLXh^h), (40) 
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where \'- k are the coefficients of lepton-up-quark-down-squark R parity breaking terms of the 
superpotential in the quark mass basis. The effective interaction (40) translates into the additional 
contributions <5Cj to Cg ; io Wilson coefficients 

V: b V ub 5C 9 = -V: b V ub 6C w = [^L\ (41) 

"QED \ m d k R J 

while no contributions are generated to C 910 Wilson coefficients [25]. For electrons in the final 
states we use bounds on A^ 2fc , X' ilk from charged current universality [59] 

( M d* \ ~ ( M d* \ 

For muons in the final state, the limits come from D + — > ir + fj, + fi~ [25]. Using the new experimental 
bound Bi(D+ -► ttW") < 8.8 x 10" 6 [10], this gives 

- - / M d« V 

The bounds on trilinear couplings (42), (43) then give the following bounds on possible enhancements 
of Cg^o Wilson coefficients for the electron or muon channel 

\V£V ub 5C e 9>w \<4A, (44a) 
\V: b V ub 6C» w \ < 17, (44b) 

with 5Cg'^ = -8Cl£. Note that in (44) the squark mass cancels. These are then added to the 
Standard Model values. The diagrams are listed on Figure 2. The possible enhancement over 
SM branching ratio predictions is quite striking and is in the case of muons in the final state by 
almost two orders of magnitude, if the values of Cg 5 io Wilson coefficients are taken to be the upper 
bounds in (44). The diagrams on Fig. 2 with photon bremsstrahlung off the final lepton pair and 
the insertion of the Q\q operator are IR divergent. We take cutoff energy to be E 7 > 50 MeV 
or Ey > 100 MeV . The contributions from various sources, the nonresonant (25) and resonant 
(28)-(30) SM contributions, the insertion of Qj, Q' 7 operator with the CV, C' 7 values given in (38a), 
and the contributions from insertion of Q9,io operators with Cq'^q bounded by (44) are summarized 
in Table 5. 

The maximal branching ratios obtainable in the framework of MSSM with R parity violation 

are 

Bl{D o ^ e + e - 7 )^ 50 MeV = 4.5 x KT 9 , Br(D° - ^V 7 )? 7 > 5 o MeV = 50 x 10~ 9 , (45) 

Br{D o ^ e + e - 7 )§ 7 > 100 MeV = 4.5 x 10~ 9 , Br^ - ^V~7)^> 100 MeV = ^ * ^ 

These are to be compared with the SM predictions (31). Note that the SM predictions are not 
affected by the cuts on the soft photon energy at the order of E^ > 100 MeV, as the bulk of 
contribution either comes from the resonances or the low p 2 region (while the cut on E^ is the cut 
on the high p 2 region). 

The enhancement due to possible R parity violating contributions is by more than an order of 
magnitude in the muon channel compared to the SM prediction. The enhancement also has a dis- 
tinct signal in the dT/dp 2 decay width distribution. In the SM model the decay D° — > either 
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Contrib. 


Br(D° -> e + e"7) 


Br(D u -► ^-7) 


Nonres. 


12.9 x lO" 11 


2.1 x lO" 11 


Reson. 


1.1 x i(r 9 


1.1 x 1(T 9 


c 7 


0.23 x 10~ 9 


0.04 x 10~ 9 


Cg 


1.37 x 1(T 9 


20.5 x 10~ 9 


Cio 


1.37 x 10~ 9 


31.3 x 10~ 9 


All 


4.52 x 1(T 9 


50.2 x 10~ 9 



Table 5: The relative sizes of various possible contributions in the context of MSSM with R parity violation. The 
photon energy cutoff is taken to be E 1 > 50 MeV. Largest possible effects are calculated. The values for nonresonant 
(Nonres.) and resonant (Reson.) LD contributions are the same as for the SM prediction. The CV denote Q7,Q'r, 
while 6*9,10 denote Qg,io insertions respectively. In the last row the maximal calculated branching ratios are given. 



proceeds through p,u>,<p vector resonances or through nonresonant two-meson exchanges, which 
are important in the low p 2 region. The R parity violating signal on the other hand would arise 
from insertion of Qgio operators and is large in the region of high p 2 (small photon energy) region 
as can be seen from Fig. 6. The largest possible effect, however, is below expected experimental 
sensitivities for rare charm decays at B-factories and CLEO-c, which are apparently expected to 
be of the order of 10~ 6 . 

4 Summary 

In this paper we have presented a detailed study of D° — ► e + e~7 and D° — > p + [i~~f decays both in 
the Standard Model (SM) and in the Minimal Super symmetric Standard Model (MSSM) with and 
without R parity violation. For the SM prediction we have carried out a calculation of the RGE 
improved Wilson coefficients of c — > u penguin operators, where Cg(m c ) has been calculated for the 
first time. The penguin operators are suppressed by V* b V u b ~ 10 -4 CKM matrix elements and are 
therefore irrelevant for the processes considered. The decays are dominated by the inclusion of Qi,2 
operators, which induce nonperturbative long distance (LD) effects. Nonresonant LD contributions 
are evaluated by employing the combined heavy quark and chiral symmetries. They are found to be 
important only in the region of low final lepton pair mass, while their contribution to the integrated 
decay width is of about 10%, or even less for the muonic channel. The decay width is dominated 
by the cascade decay D° — ► — > l + l~"f, where V = p,w, 4>. The Standard Model branching ratio 
is then predicted to be 

Br(D° - Z + r 7 ) SM = (1 - 3) x 10" 9 . (47) 

We also investigated possible enhancements of the decay widths due to new physics contributions. 
We have found that possible effects coming from gluino-squark exchanges in the context of MSSM 
with R parity conserved are masked by the LD contributions from SM. However, if the assumption 
of R parity conservation is relaxed, the tree level exchange of down squarks can increase the 
predicted branching ratios by more than an order of magnitude. The largest possible effect comes 
from the diagrams with photon bremsstrahlung off the leptons in the final state and is IR divergent. 
Choosing two different cuts on the photon energy we arrive at 

Br(D° - e+e- 7 )§ 7 > 50 MeV = 4.5 x 10" 9 , Br(D° - ^V^i^o McV = 50 x KT 9 , (48) 

Br(jD o ^ e + e - 7 )^ 100 MeV = 4.5 x 10" 9 , Br(D° - H + ^lf E ^ im MeV = 46 x 10~ 9 . (49) 



16 



Allowing for the uncertainty in the SM calculation which we discussed after Eq.(31), we consider 
that branching ratios in excess of 0.5 x 10~ 8 are not accountable by the SM. The effect of MSSM 
with R parity violation in the muon channel is the closest to the experimental sensitivities expected 
at B-factories and CLEO-c. Thus we propose the D° — > decay as a possible probe of new 

physics. 
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A Calculation of Wilson coefficients 

In this section we outline the calculation of Wilson coefficients listed in Table 1. In general, the 
Wilson coefficients at lower scale are calculated through the following steps. First the Wilson 
coefficients Ci(m\y) at weak scale are calculated by matching the effective theory with five active 
flavors q = u,d, s, c, b onto the full theory. Then the anomalous dimensions 7^) are calculated in 
the effective theory with five flavors. Using 7^, the Wilson coefficients are evolved down using 
RGE to the 6-quark scale, obtaining Ci{m b ). Then 6-quark is integrated out as an effective degree 
of freedom. This is accomplished by matching the effective theory with five flavors onto effective 
theory with four flavors. Remaining Wilson coefficients are then evolved down to charm scale using 
anomalous dimension matrices of four-flavor effective theory. Thus 



where ^4,5 (/ii, H2) are evolution matrices from scale fi2 to scale fii in four and five-flavor effective 
theories respectively, while M5 is the threshold matrix that matches the two effective theories at 
scale fx ~ rrib- 

As already discussed in section 2.1, Q10 does not mix with other operators due to chirality, so 
that Cio(// c ) = Cio(nw)- Also, the dimension five operators Qj : s do not mix into dimension six 
operators Qi,...,6 and Qg. If one is interested in these operators solely, the dimension five operators 
can be dropped from the RG analysis. We will follow this procedure and evaluate Cj separately. 
Note also, that (i) Qg operator does not mix into operators (5i,...,6 an d (ii) penguin operators Qz,...,& 
do not mix into operators Qi^- One can thus consider the RG evolution of reduced operator basis 
Qi,2, Qi,...,6 or Qi,...,6,9 if one is interested in smaller sets of Wilson coefficients C\^, Ci,...,6j or 
Ci,...,6,9 without introducing any error in the calculation. Finally, it is convenient to introduce a 
rescaled operator Qg = a/a s (uc)v-A(tt)v [34], as then the anomalous dimension depend only on 
strong coupling. 

It is instructive to do the a s counting. At leading order the RG evolution sums terms of form 
a s ln(mc//x^) which are numerically of order 0(1). At leading order one thus has to start with 
initial values Ci{mw) calculated at a®, and then evolve them using 1 loop anomalous dimensions 
(i.e. of order a s ) to get order 0(1) values Ci(fx) at lower scales. Going to higher orders an additional 
power of a s is added at each step. We thus have 



C(m c ) = U±{m c , m b )M 5 (m b )U 5 (m b , m w )C (m w ) 



(A.l) 



C i ( M ) = 0(l) + 0(a s ) + ... 



(A.2) 
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This expansion is valid also for Cg multiplying the rescaled operator Qg. Since Qg = a s /(87r)Qg, 
then Cg = 87r/a s C9, so that the expansion is 

CM = 0(l/a s ) + 0(1) + 0{a s ) + ... (A.3) 

It is thus only the NLO term that is of order 0(1) in the calculation of Cg Wilson coefficient. 
It is then consistent in a s counting to work with C9 determined at NNLO and with other Wilson 
coefficients at NLO (if one wishes to work to 0(a s )). Partial calculations at NNLO became available 
in the literature recently [42,43,60-62], however, the three-loop calculation of NNLO dimensional 
matrix has still not been performed. For this reason we will work in the following with both Cg 
and Ci e determined at NLO. 

The effective Lagrangian at weak scale n ~ my/ is 

c cS = -^L[v c * d v ud {J2 c *Qt+ E c &)+ 

* i=l,2 i=3,...,6,9 

+v c * s v us (Y,CiQt+ Yl c ^) + 

i=l,2 i=3,. ..,6,9 ^ 

+v: b v ub {Y J c t Q\+ c ^)] 

i=l,2 i=3,...,6,9 
V 1 i=l,2 i=l,2 

where Qf and the penguin and electromagnetic operators Qi, i = 3, . . . , 6, 9 are defined in Eqs. 
(3). The contributions of electromagnetic penguins have been neglected as they are suppressed 
by additional powers of a in the processes considered. In the last line of (A. 4) the unitarity of 
CKM matrix has been used V* d V ud + V* S V US + V* b V ub = 0. Above, also the masses of d, s, b quarks 
have been neglected compared to the weak scale, so that the Wilson coefficients Cj are the same 
regardless of the flavor of down quark flowing in the loop in the full theory (i.e. regardless of the 
CKM structure in front of the parenthesis in (A.4)). Thus the penguin operators do not appear in 
the effective Lagrangian at the weak scale as long as the mass of 6-quark can be neglected compared 
to m\y. This is in contrast to the case of AB = 1 decays, where up- type quarks flow in the loops in 
the full theory. Since the top quark is very heavy, its mass cannot be neglected in the loops. This 
induces penguin operators already at the weak scale. 

We start with the values of Wilson coefficients to order a s at weak scale. These are known for 
quite some time and are (in naive dimensional regularization scheme (NDR) [63]) 

Ua s (m w ) n f \ -j llajjnw) 
C 1 (m w ) = - — — , C 2 {m w ) = l-- — — , (A.5) 

while Cs^gimw) = 0. Above ^ the penguin operators do not enter the effective Lagrangian due 
to unitarity of CKM matrix. The Wilson coefficients C\^ are evolved down to fi ~ /xj using 2x2 
anomalous dimension matrix (that can be found in [63] or in Eq. (5.12) of [34]). At the scale fib the 
6-quark is integrated out, i.e. the five-flavor effective theory (A.4) is matched onto the four-flavor 
theory with the Lagrangian 

c eS = -^[v: d v ud {Y c ^+ E c i Q i ) + v c * s v us (Yc i Qi+ E c ^)] 

* i=l,2 i=3,...,6,9 i=l,2 i=3,...,6,9 , . , 

(A.6) 

= -^[v: d v ud Y,c l Qi+v: s v us Y,c l Qi-v: h v ub £ c^]. 

* i=l,2 i=l,2 i=3,...,6,9 
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The penguin operator Wilson coefficients C3 6,9 arise from the matching procedure. This is the 
only nontrivial step in the application of formulas from the literature as these were calculated for 
down-quark transitions. We use expressions for Kl — > 7r°e + e~ decay [35] where a similar matching 
procedure has to be done at charm mass, with the c-quark being integrated out. The results of Ref. 
[35] apply directly for the matching of gluonic penguin operators also in the case considered here, 
when 6-quark is integrated out, while the semileptonic Wilson coefficient C9 has to be multiplied 
by e b /e c = -1/3 • 3/2 = -1/2. We then have (see Eqs. (6.20), (8.9) of [34]) 



Zi(m b ) = Ci(m b ), 
Z 3 (m b ) ■- 



Z2(m b ) = C 2 (m b ), 
Z A (m b ) 



^F s (m b ), 



a. 



Z 5 (m b ) = ~-^ F s( m b)> 

1 a 
Z 9 (m b ) = --Z' 7V (m b ) = -lF e (m b ), 

2 47T 

with Z' 7V defined as in Eq. (8.9) of [34], while the functions 



Z & {m b ) = —F s (m b ), 



-if 



)+' 
) + ' 



Z2{p), 

(3Z 1 ( f i) + Z 2 ( f i)), 



(A.7) 
(A.8) 

(A.9) 
(A.10) 

(A.ll) 
(A.12) 



are again calculated in NDR (see also Eqs. (4.29)-(4.31) of [64]). 

The sets of operators {Qf 2 , Q3,...,6, Q9} an d {Q\ 2 , Q3,...,6, Q9} from the first line of (A. 6) are 
then evolved to the charm scale ~ m c using 7x7 anomalous dimension matrices 7 for four quark 
effective theory. The 6x6 LO and NLO submatrices involving gluonic penguins are listed in Eqs. 
(6.25), (6.26) of Ref. [34] and have been calculated in [38,65]. The remaining entries are listed in 
Eqs (8.11), (8.12) of Ref. [34] and have been calculated in [35]. 

In summary, the RG evolution from nw to fi c for AC = 1 transitions is described by the 
following procedure 



m b < fi < mw 
H = m b 
m c < ji < m b 



C(p) = U 5 (n,m w )C(m w ), 

C(m b ) -> Z(m b ), 

C{n) = Ui{n,m b )Z{m b ), 



(A.13) 
(A.14) 
(A.15) 



with U5 and U4 the 2x2 and 7x7 evolution matrices for five and four active flavors respectively. 
They can be found in Eqs. (3.93)-(3.98) of [34]. The Z(m b ) are given in (A.7)-(A.10). The values 
of calculated Wilson coefficients are listed in Table 1. 

We next turn to the calculation of effective parameters Cj\ 10 corresponding to invariant ampli- 
tudes calculated using full theory but neglecting QCD interactions. Consider c — > wj and c — ► ul + l~ 
invariant amplitudes corresponding to the diagrams on Fig. 1. The invariant amplitudes obtained 
neglecting QCD would have the same structure as one would get from the operators Q-j$ in the 
effective Lagrangian (2), if used at tree level. The parameters corresponding to these invariant 
amplitudes will be denoted C^. It is important to stress that these are not Wilson coefficients, 
as they only parametrize invariant amplitudes. They are easily obtained using calculation of Ref. 
[33] for b — > sl + l~ transitions. Following [66] we find that the coefficients are of the form 



C n = IiF I (x i ) + QiF Q {x i ), 



(A.16) 
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where I\ is connected to the weak isospin of the quarks in the loops of Fig. 1. and Q\ is their 
charge. For up quarks in the loops, as is the case for b — > sl + l~ we have // = +1, Qi = 2/3. For 
the case c — > ul + l~, that we are interested in, // = —1, Qi = —1/3, as then down quarks appear 
as intermediate states in the loops. Fj(x{) and Fq(xi) are functions of CKM matrix elements and 
masses of quarks running in the loops, Xi = m^./rriyy. The functions Fj(xi) and Fq(xi) have been 
determined by Inami and Lim [33]. Using their definitions one arrives at 



^ = -T^r- X' (A - 17) 

2 sin Vw 4 



^fsw (A18) 



with 



C = -4 *jC(xj,x d )Ii, Hi = 16 *j Sd) + 2f x d )i,] , (A.19) 

j=s,b j=s,b 

where Aj = V*jV U j/ (V* b V u b), the function C(xj,x d ) is defined in Eq. (2.14) of Ref. [33], the function 
T z (xj,Xd) in Eq. (2.7) of Ref. [33], while Fi(xj,x d ) is defined in Eq. (B.2) of Ref. [33] (note also 
the errata), the later function being changed slightly, as now 

F 1 = Q l {... }+/,... (A.20) 

(i.e. the last two lines of Eq. (B.2) in Ref. [33] are to be multiplied with /;). 

The important thing to note is that the variable Xj = m^./m^y is very small for qj = d,s,b. 
The functions C and T z are proportional to C, T z oc Xj and are thus very small. The function F\ on 
the other hand is to the leading order Fi(xj,Xd) ~ |Q/ \n{xj/xd) which is of order 0(1). Following 
the same procedure also the value of Cy L can be obtained. The leading order expressions in terms 
of Xj = nig . I rriyy are then 



Cf ~ -5/24^ \ jXj , (A.21a) 



Cf ~ -A s 16/91n (m s /m d ), (A.21b) 

C}^2^A^- i— . (A.21c) 
^ sin t/w 

The comparison of C y L coefficient and the Wilson coefficient Cg(m c ) has already been made in 
section 2.1, where it was found that Cg(m c ) ~ 10 _4 Cg L . The situation is somewhat different in the 
case of Cj Wilson coefficient. Using |U c ^U u fe| = (1.3 ±0.4) x 10~ 4 one gets \Cj L \ ~ 10~ 3 , which is an 
order of magnitude smaller than the RG improved Wilson coefficient C7. Namely, the RG evolution 
lifts the hard GIM mechanism of C7 ~ ^j x j an d replaces it with logarithmic dependence on 
the scales [i ~ m c ,mb,mw involved in the RG evolution. However, the V* b V u b suppression of Qj 
operator is still present. Again, both the inclusive rate c — > wy as well as exclusive decays are 
dominated by the inclusion of operators at one-loop level [14]. 

Finally, using Wolfenstein CKM parameters p = 0.4, i] = 0.45, A = 0.83 and the quark masses 
m d = 6 MeV, m s = 130 MeV, m b = 4.25 GeV, we arrive at = (3.9 + 1.7?) x 10" 2 . Note that (i) 
if masses of d, s, b quarks can be neglected compared to the my/, then C10 = and that then (ii) 
the low energy QCD and QED interactions cannot induce a nonzero value of C\q Wilson coefficient. 
It is thus consistent with the assumptions of OPE to set C\q = as has been done in this paper. 



20 



B List of the chiral loop integrals 



In this appendix we list definitions of the dimensionally regularized integrals needed in the evalua- 
tion of xPT and HQxPT one- loop graphs shown on Fig. 3. The integrals containing heavy quark 
propagator are 

Mm) = t£- [ d V - \ ■ = 0, (B.l) 



16vr 2 uv ' (2vr)™ J (v-q - A + id) 
j^B {0 ^ } (m,A) -J^]d Q( v . q _ A + iS){q 2_ m 2 + iS) > ( B - 2 ) 



j^ 2 C { Q^ v} (p,m 1 ,m 2 ,A) 



iff f m {l»gM»g^} < ' B ' 3 ' ) 

(2tt)" 7 - A)(^ 2 - m 2 )({q + p) 2 - m 2 ) ' 



-Y^2 D {o,^u}(Pi,P2,rn 1 ,m 2 ,m 3 ,A) = 

f d n q {h Qy, Q/iQv} 

(2tt)" y («■?- A)(g 2 - mj)((q + Pl ) 2 - m 2 )^ + p 2 ) 2 - m 2 ) ' 

where n = 4 — e. The dependence of scalar and tensor functions on v^" is not shown explicitly 
and also in Eqs. (B.3),(B.4) the iS prescription is not shown. The scalar integrals Bo(m,A), 
Co{p, mi, m2, A), Do(pi,p2, mi, mi, m 3 , A) have been calculated in [67]. We use the expressions 
of Ref. [67] in the numerical evaluation of scalar integrals Bq, Cq, Dq. The tensor integrals can 
be expressed in terms of Lorentz-covariant tensors. The notation we use for the tensor functions 
resembles closely the notation used in Ref. [68] for the Veltman-Passarino functions [69] 

B IM (m,A) = v li B 1 , (B.5) 

B^im, A) = tj^Bqq + VpVvBu, (B.6) 

Cfj,(p,mi,m 2 , A) = v^Ci + p^C 2 , (B.7) 

C flu (p,mi,m 2 ,A) = rjfa/CoQ + {v^pv + p^v v )Ci 2 + v^Cu + p^p v C 22 , (B.8) 

D /J ,(pi,P2,m 1 ,m2,m 3 , A) = v^Di + pi ll D 2 + p 2 ^D 3 , (B.9) 
D IJiU (p 1 ,p 2 ,mi,m 2 ,m 3 A) = t^Ajo + v^Du + (v^pi v + p llM v u )Di 2 

+ PlnPlvD 22 + {V^p 2 u +P2 fl V u )Di 3 +P2nP2uD 33 , 



(B.10) 



The tensor functions are calculated using the algebraic reduction [69], i.e. the tensor functions (B.5)- 
(B.10) are multiplied by four-momenta v^jp^, ... or contracted using rf" v . Then the identities such 
as v-q = v -q — A + A and/or q-p = l/2((q + p) 2 — m 2 — (q 2 — m 2 )) are used to reduce tensor 
integrals to a sum of scalar integrals. The result of this procedure has been given explicitly in [70] 
for the case of two point functions Br^^y 5 . For the case of the three and four-point functions 
C{ti,/iv}i Dip p V \ we do not write out explicitly the analytic results of algebraic reductions as the 
expressions are relatively cumbersome. For instance in the case of D^ u the final expression involves 
the inverse of a 7 x 7 matrix that corresponds to seven functions Dqq . . . D 33 appearing in the 
expression of four-point tensor function (B.10). Note as well that in this particular case there are 
ten possible relations between Dqq . . . D 33 and the scalar functions Bq, Co, Dq that one gets from 
algebraic reductions (three equations from each multiplication by v^, p% plus one relation from 

5 Note that different notation is used in Ref. [70], with B (m, A) = -7 2 (m, A)/A, B\ (m, A) = -7 2 (m, A)-Ji(m), 
Boo(m,A) = -AJi(m,A), Bu(m,A) = -AJ 2 (m,A). 
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contraction by rf*'). Obviously not all ten equations can be linearly independent. Using different 
sets of seven independent equations have to lead to the same results for Dqq . . . D33 coefficient 
functions. This fact can then be used as a very useful check in the numerical implementation. 

The aforementioned procedure runs into problems when implemented as it is in the calculation 
of D° — > / + / _ 7. Namely, for p\ = p and P2 = p + k appearing in the calculation of Cq- 4 (with p the 
four-momentum of lepton pair and k the photon momentum, see Eqs. (C.21), (C.22)) only six out 
of ten relations following from algebraic reduction are linearly independent. This problem has been 
circumvented by first calculating the tensor four-point functions with prescription k — ► k + ea, with 
a some arbitrary four-momentum and then taking the limit e — > numerically. Similarly, in the 
calculation of Cq" 5 , where p\ = k, P2 = k + p, see Eqs. (C.23), (C.24), the prescription p — > p + ea 
has been used. Because Dqq . . . -D33 are continuous functions of p\ and P2, the outlined limiting 
procedure leads to an unambiguous result. This has been also checked numerically. 

To make the paper self-contained we list in the following also the notation for the Veltman- 
Passarino functions employed by the LoopTools package [68] that has been used for their numerical 
evaluation. A general integral is 

1 N _ iff f _d^q g w • • • g^p 



16vr 2 (2ir) n J (q 2 - m 2 )({q + Pl ) 2 - m 2 ) ■ ■ ■ ({q + PiV -i) 2 - m 2 N _ i , 

with two-point functions T 2 usually denoted by letter B, the three-point functions T 3 by C and 
the four-point functions T by D. Thus e.g. Bq(p 2 , m 2 , m 2 ) and Cq(p 2 , (pi — P2) 2 ,p^,rn\,m^) are 
two-point and three-point scalar functions respectively. The decomposition of tensor integrals in 
terms of Lorentz-covariant tensors reads explicitly 

B li =p llM B 1 , (B.12) 

B^ u = V^Bqq + Pi^PivBn, (B.13) 

2 

Cp = Pllx d + p 2 pC 2 = J2Pi» C i> ( B - 14 ) 

i=l 

2 

Cp V = rjnuCoo + ^2 P'^PjvCij, (B.15) 
2 2 

Cp V p = ^ ^ {WfivPip + VvpPip + ^ppPi^CQQi + ^ ] PipJPjv p lp C i3l . (B.16) 

i=i m,z=i 
Note that the tensor-coefficient functions are totally symmetric in their indices. 



C Nonresonant LD invariant amplitudes 

In this appendix we list the analytical results for the diagrams shown on Fig. 3. They contribute 
only to the Mq U part of the invariant amplitude (12). Since separate diagrams are not gauge 
invariant, a general form of an invariant amplitude corresponding to a single diagram is 

M* = M^el(k)^u( Pl ) lu v(P2), (CI) 



<" = Ci v {p 2 )^ - C l 0kp ( P 2 )^ + Dl{p 2 )e^k aP p. (C.2) 
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For a gauge invariant sum of diagrams therefore J2i ^hq = Yli ^okp ( c -f- (13)) ^ as to ' 3e true, which 
represents a very useful numerical check. 

Note that Dq(j) 2 ) form factors corresponding to diagrams on Fig. 3 are zero. The analytical 

ji f 
y 0r/,kp v 

y1_1 



expressions for CX , (p 2 ) form factors are 



C o y = igK{v us V*B (m K , vp + A* s ) + V ud V c * d B ( m7T ,vp + A*)}, 



C, 



i_i 

Ofcp 



C, 



rn 



D 



Or] > 



(C.3) 
(C.4) 



C, 



1.2 

Or? 



-2i5ri^(K s V7 s Coo(-A;, m K ,m K , m D + A*) + KdV^Cbo(-A;, m n ,m n , m D + A*)), 



(C.5) 



/~»1_2 

u 0fcp 



/-<2_l+2_2 

^2_l+2_2 
u 0fcp 



k-p 



V...V, 



m 



us "cs 



D 



C 00 (-k,m K ,m K ,m D + A*) + (m D - vk)C 12 (-k,m K , . . . 



+ v ud v c * d 



m K -> m w , A* 



A* 



}• 



0, 



-igK 



1 



k-p 

(v-k)(v-p) m D 



■p_( 



ml- - Af )BoK, A s *) + (rn 2 , - ( mD + A:) 2 ) 

x B (m K ,m D + A^) - (m 2 K - (vk + A*) 2 ))B (m K , vk + A*) 

m K -»■ A* -> A* 



0)/ — °i 



^»2.3 



( m 2. _ (vp + A:) 2 )B (m K , vp + A s *)j + V ud V c * d 
k-p 



(C.6) 
(C.7) 

(C.8) 

(C.9) 



1 



2igK (v V* 

m 2 r, I " s cs 2(mD — vk) 



D 



Bi(m K ,vk + A*) - Bi(m K ,m D + A* 



+ J B 1 (m^,A*)- J Bi(m K ,m D -i;A; + A*) 

+ 2(m 2 K - -k 2 - A* s (vk + A* s ))Ci(-k,m K , m K , vk + A*) 

- 2(m 2 K - \k 2 - (m D - + A*)(m D + A*))Ci(-fc, m*, m*-, m D + A*) 



(CIO) 



+ v ud v; d 



2(rriD — vk) 



m K -> "v, A* -> A* J |, 



CoV = CoV( withk ^P)' 

CnV. = C V„(with k ^ p), 



^.3.2 



^»3.2 



-2U/K KsV r c * s C'oo(^, rn K ,m K , A*) + Kd^Cbo(fc, rn,,, m^, A* 

K S F C * (Coo (A;, m K ,m K , A*) - (m c - vk)C 12 (k, m K , m K , A*)) 



m 2 , 



+ KdKd( — > to,,-, A* — > A* ) 



/~<4_1 

^4_1 
u 0fcp 

/->4_2 
,^4.2 



Co 1 - 2 (with k p), 
C V p (withk^p), 
C 2 - 3 (with k <-► p), 
C 2 4(withk^p), 



(C.11) 
(C.12) 

(C.13) 

(C.14) 

(C.15) 
(C.16) 
(C.17) 
(C.18) 
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= C£- 2 (withk~p), 
= Cg(with k <- p), 



C 4 - 4 = 4i 5 ^ K s F c */„(p, k, m K , A* s ) + ^^(p, k, m n , A*) 



c okp = —ligK^ V us V*J kp (p, k, m K ,A* s ) + V ud V* d f kp {p, k, m n , A* 



k-p 



m 



°0r7 



D 

= C 4 - 4 (withk^p), 

= c Y P (withk^ P ), 

= 2igK^V us V* s - B (m K ,m D + A*) - C (p + k,m K ,m K , A*)x 
x (m^ - A* 2 ) + m D Si(m|,, m 2 K ,m 2 K ) + A*B (m 2 D , m 2 K , m 2 K ) 



+ v ud v c * d 



m K -> mm A* -> A* 



r <4.6 
^Okp 

U 0r? 



o. 



Tfl 2 Tfl 2 

^C (0,p 2 ,m 2 D ,m 2 K ,m 2 K ,m 2 K ) + —^Bo(m 2 D ,m 2 K ,m 2 K ) 



P" n ^„2 2 2 - 1 



ml, ml) + - + 



8k-p 
m*, A* A* 



2 2 

CS +5 " 2 = ^m^l^^^ColO^ 2 ^ 2 ,,^ 2 ,,^ 2 ,^!) + J?—B (m 2 D ,m 2 K ,m 2 K ) 

P 2 n /._2 2 — 2 \ i ^ 1 i T r T /"* f ___ A* A* 11 



8k-p 



B (p 2 , m 2 K ,m 2 K ) + - + KdV c * d m K — > m^, A* — > A* j j, 



y-^ V us V* s B Q (m 2 D ,m 2 K , m 2 K ) + V ud V* d B (m 2 D , w 2 , m 2 )] , 

mo* — m£>, K = ^/m£)Gjaie 3 a/(16\/27r 2 ), while in C^f 



5 3 _ iKmp 

°0»7 — = 

Cokp = 0, 



where A* = mo* — mo-, A* = 
used the abbreviation 



(C.19) 
(C.20) 

(C.21) 

(C.22) 

(C.23) 
(C.24) 

(C.25) 

(C.26) 
(C.27) 



(C.28) 



(C.29) 
(C.30) 

%t P we have 



f v (p,k,m,A) =C 0Q (k,m,m,vp + A) + (m 2 - A 2 )D 00 (p,p + k,m,m,m,A) 

- vp CqoiO 2 , k 2 , (p + k) 2 ,m 2 ,m 2 ,m 2 ) - m D C 002 (p 2 , ■■■) - AC 00 (p 2 , ...), 
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fk P (p, k,m, A) =m D Ci2{k,m,m,vp + A) + Cu(k,m,m,vp + A) 
+ (m 2 — A 2 ) Du{p,p + k, m, m, m, A) 

+ m D (D 12 (p, ...) + 2D 13 (p, ...) + D 1 (p,...)) 
+ m 2 D (D 23 (p, ...) + D 33 (p, ...)+ D 3 (p, ...)) 



m D 



-^C 00 i(p 2 , k 2 , (p + k) 2 , m 2 , m 2 , m 2 ) + -^C 002 (p 2 , • • • ) + C 222 (p 2 , ■■■) 



m 



D 

+ ^C U2 (p 2 , 
run 



m 



+ 1 + 



D 



m D J 



C 122 (p 2 ,...) + ^C 00 (p 2 ,...) 



m 



D 



+ C 22 (p 2 ,...) + — C 12 (p 2 ,...) 
m D 

- Am 2 D [C 22 (p 2 , . . . ) + Cn(p 2 , ...) + C 2 (p 2 , , 

with the dots representing the same dependence on the arguments as for the first function in the 
square brackets. 



D Nonresonant SD invariant amplitudes 

In this appendix we list the invariant amplitudes corresponding to the diagrams on Fig 2. We use 
the notation of Eq. (12), where we write down only nonzero form factors 



^ySD.I 



SD.l 



D, 



n SD.2 
SD.2 



D, 



D, 



SD.3 



m r l m 



D 



ai 



3 vk + A* 

4 k VuhV:b (fi | 1 yp c T + c! r 

3 v-k + A*\ m c J vriD a\ 
C^ D1 (withk^p), 
D^ D1 (withk^ p), 



3 vk + A* 



(*+-) 



m c run «i 



(D.l) 

(D.2) 

(D.3) 
(D.4) 

(D.5) 



SD.4 



D 



M |D.5a + M SD.56 



£ SD - 3 (with C^^Ci o 0, 
2 mo 0,1 



(') 



Z BS 



(D.6) 
(D.7) 



where A* = mo* — mn and K = ^JmjSG Fa\e 3 a /(16\/2vr 2 ) have been used, while m is the lepton 
mass. Note that the "wrong chirality" Wilson coefficients C' 7 g i0 are negligible in the SM. 
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